Renormalization Group Approach to Strong- Coupled Superconductors 
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We develop an asymptotically exact renormalization group (RG) approach that treats electron- 
electron and electron-phonon interactions on an equal footing. The approach allows an unbiased 
study of the instabilities of Fermi liquids without the assumption of a broken symmetry. We apply 
our method to the problem of strongly coupled superconductors and find the temperature T* below 
which the high-temperature Fermi liquid state becomes unstable towards Cooper pairing. We show 
that T* is the same as the critical temperature Tc obtained in Eliashberg's strong coupling theory 
starting from the low-temperature superconducting phase. A 1/A^-expansion shows that the method 
is asymptotically exact and Migdal's theorem follows as a consequence. Finally, our results lead to 
a novel way to calculate numerically, from microscopic parameters, the transition temperature of 
superconductors. 
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There is a renewed interest in understanding 
the interplay between electron-electron and electron- 
phonon interactions in strongly correlated electronic 
systems ^iiSiiiii. While the experimental data indicate 
the important role played by both electron-electron and 
electron-phonon interactions, theoretical progress has 
been limited due to the complexity associated with treat- 
ing these two interactions on an equal footing. Most ap- 
proaches rely on mean-field treatments where a "favorite" 
order parameter and broken symmetry are introduced "by 
hand". It is highly desirable if instead one could predict 
the leading instability starting from the geometry of the 
Fermi surface, the strength of the couplings, and the bare 
energy scales of the problem, such as the Fermi energy, 
Ef^ and the Debye frequency, uod- 

The renormalization group (RG) provides such an 
approach^, one that has been successful in explaining the 
stability and instabilities of Landau Fermi liquids in more 
than one dimension. Let us recall the two-stage proce- 
dure advocated there. Firstly, given a microscopic theory 
defined in all of momentum space one integrates out all 
modes except those within an energy cut-off A of Ep- In 
(i = 2, which is our focus in this work (the d = 3 case 
can be treated in analogous way^), the remaining phase 
space has the form of an annulus of radius kp (the Fermi 
momentum) and width 2A/vf (where vf is the Fermi 
velocity) . In the beginning of this process of integrating 
out high energy modes, the corrections are treated per- 
turbatively in the strength of the interactions. In the sec- 
ond stage, provided that mode elimination have reached 
an energy cut-off A <^ Efj a 1/A^ expansion emerges, 
with N c^ Ef/A. More precisely, imagine dividing the 
annulus into TV patches of size of order (A/vf)'^- The 
momentum of each fermion k is a sum of a "large" part 
(order kp) centered on a patch labeled by a patch in- 
dex i = 1, ..., A/" and a "small" momentum (order A/vf) 
within the patch^. It can then be verified that in all 
Feynman diagrams of this cut-off theory the patch index 
plays the role of a conserved isospin index exactly as in 



a theory with N fermionic species. The electron-electron 
interaction terms, written in this notation, come with a 
pre-factor of 1/N (~ A/Ef)j and the RG corrections can 
be organized in terms of powers of 1/N. Summing up 
the series of dominant corrections in order 1/N becomes 
asymptotically exact since A/" ^ oo as the RG proce- 
dure decreases the cut-off A ^ 0. This 1/A^-expansion 
has been developed in ref.[d| and has been used to show 
that, provided one can start the RG fiow at Aq <C Ef 
(large N), Landau's Fermi liquid theory emerges as the 
asymptotically exact theory for repulsive fermions with 
a generic (non-nested and with no singularities) Fermi 
surface. In this way the RG expansion goes beyond per- 
turbation theory and the end values of the couplings can 
be large. We make use of this expansion here and extend 
it for the case when phonons are present. 

Electron-phonon interactions in the Wilson-like RG of 
ref. [6], wherein momenta are rescaled to attain a fixed 
point, posed the following problem^: electronic momen- 
tum scales differently parallel and perpendicular to the 
Fermi surface while phonon momentum scales isotropi- 
cally. We circumvent this problem by using the quantum 
field theory version of the RG in which the cut-off depen- 
dence of couplings that preserves the physical quantities 
defines the fiow, with no rescaling of momenta or fre- 
quencies. 

The paper is organized as follows: The RG approach 
for interacting fermions that are also coupled to phonons 
is described in Sec. ID In Sec. [^ we illustrate our 
method for a circular Fermi surface and an analytical 
solution is obtained. The RG indicates an instability 
in the Cooper channel when the electron-phonon cou- 
pling is strong enough to overcome the effective repulsive 
electron-electron interactions. In Sec. IIIII we present the 
finite temperature formalism. We obtain the tempera- 
ture T* at which the high-temperature Fermi liquid state 
becomes unstable in the Cooper channel. We demon- 
strate that T* is the same as the superconducting tem- 
perature Tc that is obtained from the Eliashberg theory of 



strongly coupled superconductor, which approaches the 
transition from the ordered phaseifiiii. This is an alter- 
native derivation of Eliashberg's equations starting from 
a Fermi liquid state. Furthermore, by extending the 1/A^ 
analysis to the problem with phonons (Sec. II V|) . we show 
that, as in the case of the Landau's Fermi liquid theory 
for electron-electron interactions, the Eliashberg theory 
is the exact low energy effective theory obtained by using 
RG. Migdal's theorem is also derived from the 1/N anal- 
ysis. Sec. IVl contains further discussions and conclusion. 



I. RG FOR INTERACTING ELECTRONS 
COUPLED TO PHONONS 

We work in the path-integral representation and con- 
sider the general action 

5(^, 0) = 5e(^) + Sj>h{<i>) + Se-ph{^, 4>) + Se-eW (1) 

where 



Se = 






is the free electron action and 






(2) 



(3) 



is the free phonon action {k = {cc;,k} and q = {l],q}, 
where cj, Q are fermionic and bosonic Matsubara frequen- 
cies, respectively, and k, q are the momenta). The inter- 
action terms are given by 



5. 



ph= f f 9{Q)4%^ka{^,^^-,) (4) 



and 



1 T^ f I 

Se-e = :^ n / ^(^'^' ^^' ^2, ^OV^^^V^/csaV^ig ^k^a' (5) 

where /C4 = /ci + /c2 — fe. (We use units such that h = 
1 = kB') The above action defines the input physics 
at a cut-off Ao such that uod^Q^u <C Aq <C Ep (thus 
Ar-£;^/Ao>l). 

The phonons, being described by a gaussian action, can 
be integrated out exactly leading to an electron-electron 
problem with retarded interactions 

u{k4,ks,k2,ki) = u{k4,ks,k2,ki) 

- 2g{ki,k3)g{k2,k^)D{ki - ks) , (6) 



where 
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is the phonon propagator. Here we are considering 
fermions with spin and it is sufficient to focus on inter- 
actions involving two fermions with opposite spinsi. For 
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Figure 1: (a) The retarded interaction u\ (b) The self-energy 
correction; (c) The interaction vertex. 



spinless fermions the initial condition must be antisym- 
metrized with respect to 1^^2 and 3^^4. The Feynman 
diagram associated with the initial vertex u is shown 
in FigGIa). The vertex is a function of the momenta 
and frequencies of the in-coming and out-going electrons, 
with momentum and frequency conservation. The de- 
pendence on the magnitude of the momenta is irrelevant 
and all k's are taken to be on the Fermi surface. For 
a generic non-nested Fermi surface, with no singularities 
and with time-reversal symmetry, only two types of scat- 
tering exist^: forward scattering with /ci = /cs, k2 — k^ 
(these evolve into Landau parameters) and the scattering 
in the Cooper channel with k\ = — /c2, /C3 = — /C4. The 
box vertex in the diagrams in FigG]can represent either 
one of these scattering processes. The forward scatter- 
ing channel (which does not fiow under the RG as in the 
case of pure electron-electron interactions^) contributes 
to the electron self-energy i;(k, cc^) as shown in FigQIb). 
We can write 



S(cj, k) = Eo + i(\ - Z{u;, k))u; 



(8) 



with two types of contributions: a shift in the chemical 
potential (^/i ex Sq) and wave-function renormalization, 
Z(cc;,k). The shift in the chemical potential can be re- 
absorbed in the theory by assuming a fixed number of 
electrons^. The wave-function renormalization, Z(c(;,k), 
is of special interest in this problem. 

The starting point of our RG is the interaction vertex 
r*^^^[ii] in the Cooper channel, as shown in FigGJc): 



F^"^^ [u{-k3, ks, -ki, ki)] = ^t(-/c3, /cs, -ki, ki) 
u{—k, /c, — /ci, ki)u{—ks, /cs, — /c, k) 



I 



k {iuj - ek - ^{uj, k)){-iuj - ek - ^{-uj, k)) 



(9) 



where the momentum integral is such that all internal 
energies lie between and A. 



II. EXACT ANALYTICAL SOLUTION FOR A 
CIRCULAR FERMI SURFACE 

In order to establish the method and proceed with sim- 
ple analytical manipulations, we focus on an isotropic 
Fermi surface with no singular regions, and Einstein 
phonons with frequency ue (the generalization for an 
arbitrary phonon spectrum is straightforward^). Since, 
Ep^ooD^g^u we ignore radial excursions away from the 
Fermi surface in the coupling constants. In this case 
the external electron momentum can be put on kp and 
u{—ks^ks^—ki,ki) depends on the angles only via the 
difference Oi — Os. Let us define 



and 



^(^1,^3) 



-<«'/^/ 



27r 



u{-ks,ks,-ki,ki) 



where N{0) is the Fermi surface density of states. 

The RG equations are obtained from © by imposing 
the condition of cut-off independence, namely. 



d£ 



= 



(10) 



where £ = ln(Ao/A) is the RG scale. From © we obtain: 



—v{u;i,u;3,i) = 



duj Ai v{uji,uj, £)v{uj^ ^3, -^) 



,(11) 



with the initial RG condition 

v{uji,uj3,l = {)) =uq- XueD{ui -Us) 
where uq is the bare electron-electron interaction^^ and 



A = 



2N{0)g^ 

Ue 



(12) 



is the electron-phonon coupling constant. Notice that 
this is a functional equation, because retardation intro- 
duced by the phonons leads to a mixing of the couplings 
at low and high frequencies. 

The RG equation for Z^iuo)^ which for a circular Fermi 
surface does not depend on the momentum direction, can 
be likewise derived and formally integrated to give: 



Zi{u) 






Zf,uj'D{uj - oj') 



(13) 



The large- A/" considerations of ref. [d| are fully appli- 
cable here since we are dealing with non-singular interac- 
tions (albeit retarded). Thus the one-loop fiows derived 
above are exact as A/" ^ cx). More detailed analysis is 
presented in Sec. IIVI 

Notice that (fTT|) can be seen as a matrix problem: 



where 



^ = -U.M.U 



f/y(£) = v{uji,ujj,t), 



(14) 



M„-(£) = 



KiSi- 



The solution of ^^ is: 

U(^) = [1 + U(0)-P(£)]-'U(0). 
with P(£) defined as 



\i) = [ di'M{i'). 
Jo 



Therefore, the condition for the instability of the RG at 
a certain scale £ = £c is given by: 



det[l + U(0)-P(4)] =0. 



(15) 



Equivalently, we could search for the eigenvector f of 
the matrix \J~^{£c) with zero eigenvalue: 



that is. 



[i + u(o)-P(4)]-f = o, 



1 r r[u,-x..Di.-.')i^^^,^^ ^^^^ 



ttJ^J. Zl{co')co'^ + A^ 



which is an integral equation for f{uj). 

For a given value of input parameters (i^o, A, cj£;,Ao) 
the set of equations (fTH)) and (fT^ can be solved for a 
critical cut-off energy scale, Ac, at which the running 
couplings diverge and the Fermi liquid description breaks 
down. 



III. FINITE TEMPERATURE FORMALISM 

AND DERIVATION OF ELIASHBERG'S 

EQUATIONS 

The T = formalism presented in the previous section 
can be readily extended to T > 0, providing us with a 
more experimentally accessible quantity namely, a crit- 
ical temperature. We seek the temperature T* below 
which the Fermi liquid description ceases to exist as one 
scales towards the Fermi surface, that is, as Ac ^ 0. In 
this case we replace the integrals in (|T^ and (fTH)) by Mat- 
subara sums and extend the integrals in A from to 00 
to obtain 



Z{u;n)(l){Un) = -7TT*y^}uo-Xu;ED{Un- 



>(^m) 



(17) 



where we have defined (^{oOn) = f{!^n)/Z{uOn) and the 
Matsubara frequency uOn = 7rT*(2n + 1), where n is an 
integer. Moreover, from (|T^ we find: 

Z{Un) = I^XuJe y^Sgn{uJm)D{uJn - C^m) • (18) 
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Figure 2: (a) The electron-phonon vertex; (b)-(d) are self- 
energy corrections. 



The solution of (fT7)) and (fTH)) gives the value of T* as a 
function of the input parameters. 

We now relate our approach to Eliashberg's self- 
consistent mean-field theory which assumes a bro- 
ken symmetry with a superconducting order parameter 
A(a;n), in contrast to ours which starts from the Fermi 
liquid phase. Remarkably (fT7)) and (fTH)) coincide with the 
Eliashberg equations at T = Tc if we replace (piuon) by 
A(cJn)^^- This result is striking since (j){ujn) is not an or- 
der parameter and no symmetry breaking was assumed in 
our calculation. By the same token, we can show that Ac 
plays the role of the zero temperature superconducting 
gap, Aq. Since the RG procedure, approaching the insta- 
bility from high temperatures, leads to an instability of 
the Fermi liquid state at a temperature T* that is equal 
to the critical temperature Tc produced by the Eliashberg 
theory. Thus, it is no surprise that we can show that in 
the weak/intermediate coupling regime (/i* < A < 1) we 
recover the McMillan formula^ ^: 



1.13a;£;exp 



1 + A 



A-/i*(l + A) 



where 



/i 



uq, 



l^UQ\li{KQ/uJE) 



(19) 



(20) 



is the effective electron-electron interaction at the scale 
of uoe (Anderson-Morel potential) and the Allen-Dynes 
expressio] 



04. 



T* ^ {).1^^/\UJE 



(21) 



at strong coupling. 



IV. EXPANSION IN 1/iV AND DERIVATION OF 
MIGDAL'S THEOREM 

From the 1/N analysis of ref.[6], we know that the 
solution above for the one-loop RG equations contains 



the sum of all dominant corrections in 1/N . The terms 
that have been left out go to zero as A/" -^ oo. We 
have therefore demonstrated that, provided that one can 
start the RG fiow at Aq <C Ep^ Eliashberg's theory is 
the asymptotically exact description of the effective low- 
energy physics obtained by RG thanks to the small pa- 
rameter 1/N . While this is an important result, the main 
significance of our methods is that it can be used to study 
the competition between charge density wave (and other 
instabilities) and superconductivity in the strong cou- 
pling (A <C 1) regime of phonons, something that can- 
not easily be achieved in a mean-field approach. This 
issue arises for example for Fermi surfaces with regions 
of nesting. For arbitrary shapes of the Fermi surface the 
solution of the RG equations can be obtained numerically 
by discretizing the Fermi surfaces into patches^'^^. 

Since Eliashberg's theory is based on Migdal's 
theorem^^, which states that electron-phonon vertex cor- 
rections to the electron self-energy vanish as uod/Ef -^ 0, 
it must be that this theorem is built into our approach. 
We show here how this result arises and that the small 
parameter ud/Ef of Migdal's theorem is replaced here 
by 1/N. 

Let us go back to the problem before we traced the 
phonons and discuss the 1/N hierarchy. The patch in- 
dex notation carries very naturally to phonons. While 
it is obvious that when a fermion in patch i scatters to 
patch j, it emits a phonon of momentum kj — k^, we can 
also go the other way: a phonon of "large" momentum 
can be resolved, up to a two-fold ambiguity, into a dif- 
ference ki — kj associated with patches i and j, due to 
the fact that the electron momenta lie in a thin annulus 
around the Fermi surface. This fact allows us to de- 
scribe phonons with the double-index notation employed 
by t'Hooft for gluons in QCD^^: the phonon line is seen 
as made out of two counter-propagating electron lines, 
as depicted in Figj^lja). The feature of the large N ap- 
proach, that given a sum or difference of momenta, one 
can uniquely reconstruct the parts as A ^ was pointed 
out in ref.[6j. Since integrating out the phonons pro- 
duces a four-fermion interaction of size 1/N it is clear 
the electron-phonon vertex will be of size 1/\/N in our 
notation. 

Any given diagram made of n^ vertices u^ Ug ver- 
tices g and ul internal loops is of order (1/A^)^ where 
^ = — ^n — ^^/2 + TiL' This way we can organize the di- 
agrammatic expansion in powers of 1/N . The number of 
loops in each diagram, n^, can be easily obtained using 
Figl^a). Consider the problem of the RG for i;(k, cj) 
shown in Fig|2tb)-(d). It is clear from Fig. [S^b) that 
the correction of order g'^ should be taken into account. 
While g^ comes with a factor 1/A^, there is an internal 
closed loop giving an extra factor of N resulting in a cor- 
rection of 0{N^) for the self-energy. At the same order in 
perturbation expansion in 1/N ^ there are diagrams that 
are made from a simple repetition of Fig. EJb), which are 
obviously of order 0{N^) but higher order in g. Diagram 
Fig. I2Uc) is of order 0(A^^), like diagramlSlJb), and so are 




Figure 3: The numerical solution of (fT^ (continuous line) for 
II* — 0, 0.17, and 0.24 (from top to bottom) as a function of 
A. Dashed line: T*/i^E = 1.1 exp{-(l + A )/(A + /i *(l + A)} 
for /i* = 0.24; dotted line: T* /uje = 0.11 VA - 1.49 for /x* = 
0.24. The inset shows the determinant in (fT^ as a function 
of T for ji* — Q and A = 0.3, 0.4, 0.5, 0.7 and 1 (from top to 
bottom). 



all the other "rainbow" diagrams that are automatically 
included into the theory. Thus, by solving the RG equa- 
tions for a "running" self-energy at one-loop, we in fact 
take into account all corrections of order 0{N^) to all 
loops. The infinite series of diagrams being summed is 
not arbitrary and arises naturally from the RG equations, 
as shown by our 1/A^ analysis. The final diagram of order 
g^ is shown in FiglS^d): this is the famous vertex correc- 
tion to the electron self-energy due to electron-phonon 
interactions studied by Migdal^^. Notice that diagram 
EJd) does not contain any internal loops and is of order 
0{N~'^) and thus is vanishingly small as A^ ^ oo. It 
is easy to show that these corrections to the self-energy 
lead to (|T3)) . Thus, Migdal's theorem is built into the 
large N approach^^. We can also show that the phonon 
self-energy has only contributions of order 1/N'^. 



V. DISCUSSION AND CONCLUSION 

Our method provides a novel way to calculate Tc of su- 
perconductors, starting from the normal state. At the su- 
perconducting instability, Eliashberg's equations are de- 
rived from the RG fiow equations. These self-consistent 
integral equations (|17|) are equivalent to the condition 
of zero determinant expressed in (|T^ . In our RG ap- 
proach, once the bare values of the couplings are given, 
we simply calculate the fiow of det[l + U(0) • P(T)] by 
solving the differential RG fiow equations and determine 
T* as the point when this quantity becomes zero (see in- 
set in FiglSJ. This numerical calculation is much simpler 
than solving (|T7)) directly, since it does not involve self- 
consistency conditions. In FigElwe show the calculation 
of T* as a function of A with this new method together 
with the approximate asymptotic expressions for weak 
and strong coupling. 

In summary, we have developed an asymptotically ex- 
act RG method that takes into account electron-electron 
and electron-phonon interactions in an unbiased way and 
reproduces all of Eliashberg's theory and also provides 
a framework (large- A/") for understanding it as well as 
Migdal's theorem that goes into it. Our procedure can 
be used for any Fermi surface geometry and for any num- 
ber of scattering channels (forward, charge and spin den- 
sity wave, etc) and therefore allows for the study of the 
competition between scattering channels^. Finally, our 
procedure allows for a new numerical way to investigate 
superconductivity in metals. 
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